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Abstract 

The dynamics of a binary system with non conserved order parameter under 
a plain shear flow with rate 7 is solved analytically in the large- iV limit. A 
phase transition is observed at a critical temperature T c ("y). After a quench 
from a high temperature equilibrium state to a lower temperature T a non- 
equilibrium stationary state is entered when T > T c (j), while aging dynamics 
characterizes the phases with T < T c {p/). Two-time quantities are computed 
and the off-equilibrium generalization of the fluctuation-dissipation theorem 
is provided. 

05.70.Ln; 64.75.+g; 83.50.Ax 
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I. INTRODUCTION 



In the field of modern statistical mechanics many efforts have been devoted to understand 
the properties of non-equilibrium. These studies were promoted by the recognition that 
many relevant systems, among which glasses and coarsening systems, are intrinsecally far 
from equilibrium and, therefore, cannot be described by theories formulated as extensions 
of Gibbs statistical mechanics. Instead, specific concepts relative to the non-equilibrium 
state must be developed. A promising subject in this field is the attempt to build an 
out of equilibrium fluctuation-dissipation theorem (FDT). Specifically, the aim of these 
studies is to investigate if the relation between the (integrated) response function \ an d the 
autocorrelation function D is meaningful and bears some relevant informations regarding 
the non-equilibrium state. There are several reasons for expecting this, the most natural 
being induction from the equilibrium case. In equilibrium these two quantities are linearly 
related by means of the standard FDT and a fundamental quantity, the temperature of the 
system, appears as the coefficient of this relation. For off-equilibrium systems we do not 
have nowadays an analogous theorem, of such a wide generality, but conjectures, indications 
or, sometimes, proofs exist claiming that, in restricted contexts, the relation between x an d 
D can be inferred on general grounds. In particular, the slope of the curve x(D) has been 
interpreted as an effective temperature of the system, different from that of the thermal 
bath. 

These concepts, which were originally developed in the study of turbulence [1] have been 
successively applied to the glassy state [2,3] and, more generally, to aging systems [4] where 
the non-equilibrium state is generally determined by the change of a thermodynamic control 
parameter, as in a temperature quench. 

A different non-equilibrium situation is very often considered because of its physical 
relevance: The case of driven systems. These systems are subjected to an external forcing 
that prevents them to attain equilibrium. Examples of this situation range from systems 
subjected to a force which does not derive from a potential, as sheared fluids or conductor 
wires with a potential difference applied at their hands, to systems with externally imposed 
thermal gradients or magnetic materials under the action of an oscillatory magnetic field 
and tapped materials. Despite many studies on this subject [5] much less is known on the 
effect of the drive on the fluctuation-dissipation relation, despite some progress has been 
done recently [3,6]. 

The complexity of the glassy state, and the difficulty of both analytical and numerical 
results in this field, make the task of understanding the mechanisms whereby the FDT can be 
generalized a hard question. Therefore, in this paper we focus on binary coarsening systems 
which, besides being physically interesting of their own, can be regarded as simple paradigms 
of aging behaviour. We solve exactly the dynamics of a model for a ferromagnet quenched 
from an initial disordered state under the action of an external driving field. This kind of 
process is relevant in many physical contexts, as it will be discussed later. We consider a 
coarse-grained system described by the usual double-well ip A hamiltonian where is a scalar 
order parameter field. A Langevin dynamics is considered and the order parameter is not 
conserved. The driving field is a planar shear or Couette flow. This model, without the 
shear flow, is referred to in the literature as model A. The choice of this particular flow, 
besides analytical tractability, is also motivated by its practical relevance, which appears 
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clearly in the physics of complex fluids and binary mixtures, where the behaviour of the 
stress as a function of the shear rate is one of the main characterizations of fluid systems. 
Moreover, the structure of binary and complex fluids is strongly affected by the flow and 
this is relevant in many applications [7] and interesting for its theoretical implications. In 
phase separation of binary mixtures, for instance, it was found that not only the morphology 
but also the laws of dynamical scaling are affected by the shear-induced anisotropy [8,9] and 
new phenomena related to the presence of stress in segregating systems, as logarithmic-time 
periodic oscillations [10] and violation of dynamical scaling [11] have been observed. 

In this paper we focus on a non-conserved order parameter; however we believe that our 
study could be also useful as a preliminary step for understanding more complex systems 
with conserved ordered parameter. By its own the model of this paper can describe the 
behaviour of liquid crystals undergoing the isotropic-nematic transition where the order 
parameter is not conserved [12]. The phase separation properties of this model have been 
studied in [13] where only the behaviour of one-time correlations was considered. 

The role of the non-linear term in the y? 4 hamiltonian is crucial for quenching below and 
at the critical temperature. However its presence makes the analytical study of the model 
impossible unless one considers approximate theories. Among these the so-called large- 
N limit has played a prominent role in the study of phase separation. The approximation 
takes into account a vectorial system with an infinite number of components N. This model, 
which can be solved by a self-consistent closure, is one of the few cases where it is possible 
to illustrate the out-of-equilibrium behaviour of a coarsening system by exact calculations 



In the case without shear the large- N approximation captures the essence of the phe- 
nomenon at a semi- quantitative level [14]. In this work we use this approximation to calculate 
the two-time correlation function showing how the dynamical exponents are affected by the 
presence of the flow. The response function can be also exactly computed allowing a careful 
discussion of the behaviour of the fluctuation-dissipation relation. 

The plan of the paper is the following. In Sec. II the model is defined and in Sec. Ill 
its stationary critical properties are discussed. In Sec. IV the self-consistent condition is 
explicitly worked out and one-time quantities are computed. Two-time quantities, the au- 
tocorrelation function and the response function, are studied in Sec. V and Sec. VI. The 
violation of the FDT is considered in Sec. VII and a final discussion is presented in Sec. VIII. 
Appendixes A-D, containing some mathematical details, complete the paper. 



where {f a } are the N components of the vectorial order parameter. For instance, in the 
case of magnetic systems, <p is the local magnetization, T is a transport coefficient and the 
gaussian white noise 77, representing thermal fluctuations, has expectations 



[12]. 



II. THE MODEL 



We consider a system with dynamics described by the equations 






(r) a (x,t)) = 

(r) a (x, t)r)P(x', t')) = 2T5 a pT5(x - x')5(t - t') 



(2) 
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where T is the temperature of the thermal bath. Without the convective term on the 
left hand side, Eq. (1) would be the usual Langevin equation which governs the purely 
relaxational dynamics of a system with non-conserved order parameter and hamiltonian 
H[<f\. In that case the relations (2) would assure that in thermodynamic equilibrium at 
temperature T the fluctuation-dissipation theorem is verified. Here Eqs. (2) are supposed 
to hold on the basis of local equilibrium [5,15]. More precisely it is assumed that the 
non-equilibrium system can be subdivided into cells small enough that any thermodynamic 
property - which in nonequilibrium situations may depend on space - vary slightly over 
one cell, but large enough that equilibrium statistical mechanics holds within each cell [16]. 
Furthermore, since thermodynamic properties out of equilibrium may depend on time, the 
characteristic time over which a macroscopic fluctuation dies away within one cell must be 
much smaller than the typical evolution time of the system. This implies that over intervals 
much smaller than this evolution time local equilibrium is mantained in each cell and the 
fluctuation-dissipation theorem holds. This will be recovered in Sec. VII. 
The velocity v in Eq. (1) is chosen to be the steady planar shear flow 

v = >yye x (3) 

where 7 is the spatially homogeneous shear rate and e x is a unit vector in the flow direction. 
We observe that, in spite of the presence of a velocity field proportional to a space coordinate, 
translational invariance still holds, since a shift of a in the y-direction is equivalent to a 
galilean transformation to a new reference frame globally moving with an added velocity 
07 in the x-direction. This allows in the following to introduce Fourier transforms and the 
standard definition of the structure factor [8,15]. The hamiltonian has the Ginzburg-Landau 
form 



H[<p\ = I d d x 
J v 



1 I V y?| 2 +V + — (y? 2 ) 2 

2 I Sk 4tv v ^ ; . 



(4) 



where r < 0, g > and V is the volume of the system. We will be interested in processes 
where the system is initially in the equilibrium (7 = 0) infinite temperature state with 
expectations 



< ^ a (x,0) >= 

< <p a (x, 0)</(£', 0) >= A 5 a/3 5(x - £') 



and then is suddenly put in contact with a heat bath at temperature T, and subjected to 
the shear flow. A is a constant. 

In the large-N limit [12] the equation of motion for the Fourier transform <p a (k,t) = 
j v d d xip a (x, t) exp(ik ■ x) takes the linear form 

" 7fe^S^ = ~nk 2 + I(t)]v a (k,t) + V "(k,t) (6) 



dt dk y 



where 



(V a (k,t)} = (7) 
( V a (k,t)^(k',t')) = 2T5 a pVV5^5{t - 
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and the function 



I(t)=r + jj<<?(x,t)> (8) 



has to be calculated self-consistently. The quantity < (^{x^t) > can be expressed as 

±<0\x,t)>=±Y.C(k,t) (9) 

k 



where the structure factor 



C$,t) = j^<<ftk,t)>0(-k,t)> (10) 
is solution of the equation 

- 7**^j|^ = -2T[e + I(t)]C(k, t) + 2YT. (11) 
The momentum sum in Eq. (9) extends up to a phenomenological ultraviolet cut-off A. 



III. STATIONARY PROPERTIES 

Letting dC(k,t)/dt = in Eq. (11), and setting T = 1, the stationary form of the 
structure factor reads [15] 

C(k) = 2T H e-'fo^'^'-^dz , (12) 
Jo 

where 

IC(s) = k + ^sk x e y (13) 

and 

C 2 = lim/(f)=r+ (14) 

k 

— * 

This quantity plays the role of a transverse correlation length, because the modes of C{k) 
with k x = take the usual Ornstein-Zernike form 



T 

Inserting Eq. (12) into Eq. (14) one has 



C(k x = 0,k ± )= J2 2 . (15) 



Zf = r + ^2T^2 J°° e- 2 ^ K2{t ' )dt '- 2z ^dz . (16) 

k 

— * 

This equation can be conveniently solved separating the k = term under the sum. Then, 
for very large volumes 
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C 2 = r + gTP{^) + g T (17) 

v s± 

where 

P(C 2 ) = 2 / 7^e- fc2 /A 2 f ^/W*-*^ . (18) 
J {27r) a Jo 

The existence of a microscopic lengthscale proportional to A^ 1 corresponds to the presence 
of a microscopic relaxation time 

r M = (2A 2 )" 1 (19) 

which has to be compared with the shear flow timescale 

T S = T 1 . (20) 

As usually, it is assumed [8,16] that the flow does not distort the structure of the system at 
microscopic level so that tm <C t s and 

2A 2 

A = > 1 . (21) 

7 

In the following we will concentrate on the particular cases d = 2 and d = 3. The function 
P(x) is a non negative monotonically decreasing function with the maximum value at x — 0. 
It is related through Eq. (70) to the function f L (x) calculated in Appendix A and, from 
Eqs. (78,85), its value at x = is given by 

where T(3/4) is the Gamma function evaluated at 3/4. By graphical analysis one can easily 
show that Eq.(17) admits a solution with a finite value of £j_ for all T. However, there 
exists a critical value of the temperature T c {^) defined by 

r + £T c ( 7 )P(0) = (23) 

such that for T > T c (7) the solution is independent of the volume, while for T < T c (7) it 
does depend on V . From Eqs. (22,23) one has 

f (47rM 2 )/ln(A) , d = 2 



Tc(7) \T c ( 7 = 0)[l-r 2 (3/4)3- 1 / 4 v /2T7A] _1 , d = 3 (24) 

with Ml = -r/g and T c (7 = 0) = 4M 2 7r 3 / 2 /A. Notice that the effect of the flow is to 
increase the value of the critical temperature with respect to the case with 7 = [17], 
according to T c {^) ~ In 7 in d — 2 and Tj^j) — T c (0) oc in d = 3. In d = 2 this produces 
a qualitative difference because a phase-transition occurs at a finite temperature, differently 
from the unsheared case. 
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Finally, concerning the behaviour of the transverse correlation length, above T c (pj) with 
[Tp — T c ('j)]/T c ('y) <C 1, from Eqs. (17,18) and the expression of f L calculated in Appendix 
A (Eqs. (77,84)), one finds 

C 2 = M 2 ( ^M ) (TA 3 + 1/g)- 1 d = 3 (25) 

where A 3 and B 2 are defined in Appendix A. Hence, for the exponent v defined by £± ~ 
( P ^(i) > )~ 1 ' one fi n ds ^ = 1/2, with logarithmic corrections for d = 2. At T c (7), one has 
£I 2 ~ 1/v 7 ^ 7 in d = 3 and £j 2 | log(^J 2 )| ~ 1/V for d = 2. The behaviour of £i for T > T c ( 7 ) 
is shown in Fig.l. 

For fixed 7, £j_ decreases by increasing the temperature. Indeed this is expected because 
coherence is suppressed by thermal fluctuations. The role of 7 is more subtle. Naively one 
would expect that increasing 7 produces a smaller because coherent island are washed 
away by the flow. However it must be also considered that 7 raises the critical temperature 
T c ( 7 ). Hence, for fixed T, larger values of 7 make [T — T c ( 7 )]/T c ( 7 ) smaller, the system gets 
closer to the critical point and the coherence length is increased. This second effect competes 
with the former and produces a net increase of £j_ with 7. However, by changing 7 and T in 
order to mantain a constant distance from criticality, namely keeping [T — T c (j)]/T c (^) fixed, 
the second effect is suppressed and one sees that, indeed, £j_ is reduced by the flow (see inset 
of Fig. 1). For small [T — T c ('j)]/T c ('j) the power law behaviour (25) is observed. Finally, 
for large temperatures, £j_ becomes 7-independent, because the de-coherence produced by 
thermal fluctuations prevails over shear effects. 

Finally, turning to the phase T < 2^(7), the transverse correlation length is given by 
ei 2 = T/\m 2 V) where 



M=M T~m (27) 



is the analogue of the spontaneous magnetization in equilibrium. 



IV. DYNAMICAL SOLUTION 

We will solve the model at asymptotic times in the case of quenching processes with 
initial conditions given by (5). From now on we will take the infinite volume limit. The 
formal solution of Eq. (11) is given by 

C(k,t) = 4re- 2/ ° duK2 ^e- k2 ^ 2 + ^ f e -* e -# I* y{z)dz (28) 

1 ; y(t) y(t) Jo 

where 
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FIG. 1. is shown against T for different values of 7 in d = 3. The inset shows the same 
function plotted against [T — T c (~f)\/T c {j). The reference line below the data curves represents the 
expected slope 1/2 (Eq. (25)). 



e- k2 ^C(k,t) 



y{t)=exp (2^ [r + gS(t')]dt 

and, from Eq. (9) 

S « = fe^fWW = & ^CAt) = / g 
Then, denning /(t) as 

f (f) = /* c -2 / '[/c 2 W]d*-fc 2 /A 
1/17 7 (27r) d 



+ 1/2^)^2 (4 - + !_^) 



(29) 



(30) 



(31) 



(32) 



one has 
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S(t) = -L ( Ao /(t) + 2T jf* /(* - f )l/(f • (33) 

We observe that the correction induced by the flow in Eq. (32) is independent on the dimen- 
sionality of the system and that, in absence of flow , the expected behaviour f(t) (8irt)~ d / 2 
is recovered. 

From Eqs. (29,33), one obtains the integro-differential equation 



dy(t) 
dt 



2y(t)(r + gS{t)) = 2y{t)r + 2g (Aq f(t) + 2T J* f(t - t')y(t')dt^ (34) 



that can be solved introducing Laplace transforms denoted by f L (s) = J °° f(t)e~ st dt and 
assuming that y L (s), the Laplace transform of y(t), is well defined, as it will be verified a 
posteriori. 

Eq. (34) implies that y L (s) can be expressed in terms of f L (s) through the equation 

v L (s) - l + 2A ^ (35) 

The function f L (s) has been calculated in Appendix A. Then, from Eq. (35), inverting the 
Laplace transform (see Appendix B for details), the function y(t) can be obtained. For 
t >> r s , y(t) behaves as 



d=2 



y(t) 



1 (1 , M o A o \ 
2B 2 T\2 ^ 2T C ^)> 

1 ( 1 



+ 



log(2£ 



M l£o\_ 1 



-) 6 



2B 2 T C \2 1 2T c ( 7 )/logt 

H(T + A M 2 )r 2 



T > T c ( 7 ) 

T = T c ( 7 ) 
T < T c ( 7 ) 



(36) 



T > T c ( 7 ) 

T = T c ( 7 ) (37) 
T < T c ( 7 ) 

where M is given in Eq. (27), £j_ and T c ( 7 ) are given in Eqs. (24,25,26) and the parameters 
A 3 , B 2 , B 3 are defined in Appendix A. Eqs. (36,37) hold for t > in the high temperature 
phase T > T c ( 7 ) and for t > r s when T < T c ( 7 ). 

V. AUTOCORRELATION FUNCTION 

In this Section we analyze the asymptotic behaviour of the autocorrelation function 



d=3 



y(t) 



l r -2£7 2 t l+MgAo/r e ( 7 ) 
2 C 1+4T C ( 7 )A 3 

1 l+M 2 A /r e 

2 1+4T C ( 7 )A 3 
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where the correlator 

V(k, t; k', t') =< ip(k, t)ip(k', t') > (39) 

satisfies the equation 

dV(k,t;k',t') , dV(k,t;k',t') ri2 „, , -.. 

1 ' j ' ~ iK 1 ' ' ; = —[k + gS(t) + r}V(k, t; k , f) 40 

at oky 

with initial condition V(k, t'; k',t') = C(k' ,t')5(k + k') and * > t'. 

The function D(t,t') is calculated in Appendix D using Eqs. (36,37) for y(t). Here we 
discuss the results in the three cases T > T c (7), T = T c ( 7 ), and T < T c (t)- 



T>T C ( 7 ) 

In the regime £' > ^ Eqs. (115,117) show that the correlation function becomes the TTI 
quantity 

AT fOO _n , 1 

(41) 

with t — t — t'. This implies that £ 2 is the characteristic relaxation time of fluctuations 
above the critical temperature. 



T = T c ( 7 ) 

At the critical temperature, due to Eqs. (118,120), D(t,tf) tends, for t' > t s (Eq. (20)), 
to the form D st (r, oo), namely 



4-/ r°° , ■ 

D(t,tf) ~ D„(T,oo) = — m (y + l/2A'T"r- dy. (42) 



Then one has TTI again but, for large r, D st (r, oo) decays as a power law D st (r, oo) oc 
r~ rf / 2 , implying the absence of a characteristic relaxation time, as usual at criticality. Notice 
that the shear has the effect of increasing the exponent of the power law decay with respect 
to the value — (d — 2)/2 of the undriven case. In other terms, the flow tends to decorrelate 
faster the system. 



T < T c ( 7 ) 
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In this case there are two time regimes of interest: 

i) short time separation: t' — > oo , ^7^0 (Quasi-stationary regime) 

ii) large time separation: t' — > oo , ^ — > oo (Aging regime). 
In the time sector i), from Eqs.(124) one has 

D(t,t') = M 2 + D st (r, oo) (43) 

where D st (r, oo) is the same time translational invariant quantity found at T = 2^(7), 
Eq. (42). 

On the other hand, in the limit ii), from Eq. (125), one gets 

( d+2 d+2 
7) 4 (l + t/t)-*? . 2 ^ (44) 

Y i+t'/t (l+t'/*) 2 

Notice the dependence of D ag on the ratio t/t' alone, as usual in slowly relaxing aging 
systems. Furthermore, for t/t' 3> 1, a generalization of the Fisher-Huse exponent A defined 
by 

D(t,t') oc (t'/t) x (45) 

gives 

A = 1±». (46) 

The structure of the autocorrelation function provided in Eqs. (44,43) for large t' allows 
its splitting into the sum of two distinct contributions 

D(t,t') = D st (r,oo) + D ag (t/t'). (47) 

In the regime i) D ag {t/t') remains fixed to its initial value D ag (l) = M 2 , while D st {r, 00) 
decays to zero. This makes Eq. (47) consistent with Eq. (43). Conversely, in the regime ii), 
D st (r, 00) has already decayed to zero and the whole time dependence is carried by D ag (t/t'), 
providing again consistence between Eqs. (47) and (44). 

This whole pattern of behaviours of D(t,t') is closely reminiscent of what is known for 
the system without shear [18]. In that case it was shown that this features reflect an un- 
derlying structure of the order parameter field which, in the late stage of the dynamics of 
a quench below T c , can be decomposed into two statistically indipendent stochastic com- 
ponents, (fi(x,t) = ip(x,t) + a(x, t), responsible for the stationary and aging part of the 
autocorrelation function respectively, as discussed in [19]. In the present case, by proceed- 
ing along the lines of the case 7 = 0, it is straightforward to show that the same splitting 
of the order parameter field can be explicitly exhibited, which accounts for the structure of 
D(t,t') discussed above; we refer to [18,20] for further details. In systems with a scalar or- 
der parameter [21] the physical interpretation of this decomposition is the following: given a 
configuration at a time well inside the late stage scaling regime, one can distinguish degrees 
of freedom ip(x,t) in the bulk of domains from those a(x, t) pertaining to the interfaces. 
The first ones are driven by thermal fluctuations and behave locally as in equilibrium; the 
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second ones retain memory of the noisy initial condition, and produce the aging behaviour. 
For large- N, where topological defects are unstable and domains are not well defined, the 
recognition of the physical degrees of freedom associated to ip(x, t) and a(x, t) is not straight- 
forward as for scalar systems. Nevertheless, the possibility of splitting of the order parameter 
also in this case indicates that the same fundamental property is shared by systems with 
different N. 

The overall behaviour of the autocorrelation function is shown in figure 2, for different 
t' . This figure shows also a comparison with the case 7 = 0. The pattern is qualitatively 
similar in both cases: Initially a fast decay to a plateau value is observed. This is due to 
the vanishing of D st (r, oo) in the quasi- stationary regime while D ag (t/t') remains constant 
D ag (t/t') ~ M 2 (this value corresponds to the height of the plateau in the figure). Notice 
that curves with different t' collapse in this regime, due to the stationarity of D st (r, oo). For 
longer times, when t — t' ~ t', the autocorrelation function departs from the plateau value, 
because also D ag (t/t') starts to fall off, and a pronounced dependence on the waiting time 
t' , or aging, is observed. After the stationary regime is over (around r ~ 10) also 7 plays an 
important role. Actually one observes that the decay of D ag (t/t') is promoted by the shear, 
as expected because correlations are washed out by the flow. Indeed, considering the case 
with t' = 10 5 , for example, D(t,t') decays from one to the value 10~ 3 in a time r of order 
10 9 for 7 = and in a time of order 3 • 10 7 for 7 = 0.1. For long times this effect results in 
the larger Fisher-Huse exponent (46) with respect to the undriven case, where A = d/A [12]. 



VI. LINEAR RESPONSE 



In order to calculate the response of the field to an external perturbation h(x, t) we add 
to the original hamiltonian the term 5Ti — — J dx d h(x,t) ■ <p(x,t). The Langevin equation 
(1) for a generic component becomes 



d(p a (x, t) 



+ V-( V a (x,t)v) 



5H[0\ 



dt ^ v ' ' ' 5<p a (x,t) 

The two-time response function is defined as 



+ h a (x,t)+T] a (x,t). 



(48) 



TZ(k,t;k',t') = (2tt) 



d 6<<p a (k,t)> 



5h a (-k',t') 



(49) 



{h(k',t')=0} 



with t > t' . Here, due to rotational symmetry, a generic component of the order parameter 
can be considered and vectorial indices can be dropped. Using the properties < (f a (k, 0) >= 
and < i] a (k, t) >= we get 



K(k,t; k',t') = (2ir) d 5 (£(t - t') + k') , ^e-f duK ^ 

v ' \ y{t) 



(50) 



where the function y(t) is defined in Eq. (29). The autoresponse function can be also 
introduced as 
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FIG. 2. Autocorrelation function vs r for different waiting times t', with T = T c ( / y)/2, d = 3, 
A = 1. The upper figure refers to the case without shear, the lower one to the case with 7 = 0.1. 
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\ y(t) J 

2 l - d ' 2 I y(t') ft - t' 
(4tt)^ y(t) { 2 



)/A 2 



-d/2 



+ r M 



4 + 4 2 (t-t>)3 _ 2 ft~^) 4 
^ ^ 3 ' t-t'+2r M ' (t-t'+2r M ) 2 



(51) 



where = /C 2 (t - f) - k 2 . 

Computing the integrated response one obtains 



X (t,t w ) = f R(t,t')dt' = 

Jt w 



J dzy/y(t-z)- 



-d/2 



4 + h 2 z 2 



(52) 



where we have introduced the variables z = t — t', r = t — t w and neglected the two 
asymptotically irrelevant terms containing tm in the square root in the last line of Eq. (51). 
We will evaluate this function in the asymptotic regime t w — > oo using the large time 
behaviour of y(t) given in Eqs. (36,37). 



T > T c ( 7 ) 

Inserting Eqs. (36,37) in Eq. (52), the response function reads 



Xst(r, fj 



(4vr) rf / 2 



Jo 



2z/2 (z + 2r M )^ 

2 -,2 



(53) 



which is a time translational invariant quantity, as expected, because the system reaches a 
stationary state. 



T = T c ( 7 ) 

For t w 3> t s one has 



X(t,t w ) = Xst(r, oo) 



dz 



[z + 2r M ) 



-d/2 



(Any/ 2 Jo ^1 



zl 2 z 2 



(54) 



showing that also at criticality the response function is time translational invariant. 

After an integration by parts, the function Xstiji oo) can be expressed in d = 3 in terms 
of the 2 -Pi hypergeometric function as 



Xst(r, oo) 



1 



(2tt) 3 / 2 
9 



A (l/2 + (A 2 r + l)- 1/2 (4 + I/37V)- 1 



T 6 ' 2 2 F 1 l/2,3/4,7/4,- 7 V/12 



/2 



(55) 



while, in d — 2, one gets 
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Xst(r, oo) 



1/(2tt) 



\A + V(3^ 2 ) 



log(A 2 r + 1) + log 1 + Jl + 1/(3A 2 ) 



- log (1 + ^/1 + 1/(3A 2 )^1 + l/3( 7 r/2) 2 - yr/(6A) 
In particular, for r <r s , one has the limiting behaviour 



(56) 



^log(A 2 r + l) d = 2 
_JL(l i_\ d - 3 

r c (7) v 1 v^nJ a ~ d 



X st (r,oo)= { ^ , , , , (57) 

while 



= (58) 



in the opposite limit r ^> t s 



T < T c ( 7 ) 

For T < T c ( 7 ) the integrated response function is given by 

y(tt) = [ T ^ u + 1 ^y d/2 {] - u/tYW (59) 

The term (1 — u/t)~( d+2 ^ 4 can be expanded in powers of w/t and integrated term by term. 
It can be shown that the first contribution is dominant for t w ^> t s and behaves as Xstij, oo) 
in Eq. (54). Then in the limit t w 3> r s the response function becomes TTI also in the low 
temperature phase. 

It is interesting to study the behaviour of the first correction to this result, namely the 
next term appairing in the power series in u/t discussed above. The dependence on the two 
times of this quantity, that will be denoted as Xag(t,t w ), can be factorized as Xag(t,t w ) = 
t _1 6(r). Here b(r) is a TTI function that, starting from zero at r = 0, saturates to a constant 
value for r > t s . In the regime r > r s , therefore, one has Xag{t,t w ) ~ In conclusion, 
taking into account also this first correction, it results that x(Muj) — Xstiji oo) + Xa g (t,t w ) 
and, recalling the behaviour of Xag(t,t w ), one concludes that for large t the second term is 
always negligible with respect to the first, as anticipated after Eq. (59) 

A similar structure, with a stationary and an aging part, is also found in the case with 
7 = 0. However, in that case, one finds a similar result, namely Xag(t,t w ) = t _1 6(r) and the 
properties of b(r) discussed above, only above the upper critical dimension du = 4. Instead, 
below du it is Xag(t,t w ) = t~ a B(t/t w ), for r/t w > 1, with a = (d — 2)/2. Notice that this 
form implies that Xagit, t w ) does not vanish for large t at the lower critical dimension d = 2. 
A similar pattern is also found for scalar systems [21]. 

When shear is applied the system behaves as in the case with 7 = and d> djj. This 
is due to the fact that the mathematical structure of the solution of the model with applied 
flow in dimension d is similar to the structure of the solution without shear in an increased 
dimensionality d + 2. This can be checked, for instance, from the behaviour of the function 
y(t) of Eqs. (36,37) as compared to the form without shear [18]. This phenomenon raises the 
effective dimensionality at or above du = A for the cases d — 2,3 considered in this paper. 
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VII. FLUCTUATION-DISSIPATION RELATION 

In this Section we will discuss the out of equilibrium generalization of the FDT, namely 
the relation between D(t,t w ) and x(t,tw) i* 1 the large t w regime. 

T > T c ( 7 ) 

Let us first recall that in the case without shear, for times t w larger than the equilibration 
time, D(t, t w ) and x(t, t w ) become TTI quantities D{t) and x( r )- Moreover, since the system 
attains an equilibrium state, standard FDT holds 

T X (r) = D(r = 0)-D(r). (60) 

The effect of the shear flow is to transfer energy to the system, and, also in the quench 
to a temperature above the critical point where TTI is obeyed, the equilibrium Gibbs state 
is never reached. Then the FDT does not hold and D(t,t w ) and do not satisfy an 

equation analogous to Eq. (60). The fact that D(t,t w ) = D st (r,^±) and x(t^w) = Xst( r , £±) 
are TTI quantities implies that the dependence on the two times of x{t>tw) can still be 
accounted for through D(t,t w ), yielding a non trivial relation x(D) between an d 
D(t,t w ). The fluctuation-dissipation plot, namely the relation x(D), is shown in Fig. 3. 

1 
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FIG. 3. Tx(t) is plotted against D(t) for a quench to a final temperature T = (21/20)T c ( , y), 
with t w = 10 5 , d = 3 and A = 1. Solid lines correspond to 7 = 10~ 2 , 10 -1 , 1, from top to bottom. 
The dashed line, which is almost indistinguishable from the case 7 = 10 -2 , is the equilibrium 
behaviour 7 = 0. The inset shows the fluctuation-dissipation ratio X(D) defined in Eq. (62). 

The plot follows the straight line (60) for small values of r, when r ^ r s . Recalling 
the discussion of Sec. II regarding local equilibrium one has that the system behaves as in 
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equilibrium up to timescales of the order of the evolution time of the system which, in this 
case, is the flow timescale t s . 

For larger values of r (smaller D) x{D) deviates from the equilibrium form (60) due to 
the effects of shear. In particular, for T > ^(7) its asymptotic (r = 00) value is Tx(r, ~ 
M 2 T/T C (2 7 ) (Eq. (58)) which is less than M 2 — M 2 = M 2 T/T C ( 7 ), the value obtained from 
Eq. (60). From Eq. (58) it is readily seen that the relation between the response x{t, t w ; 27) of 
a system subjected to a shear flow with rate 27, and the autocorrelation function D(t, t w ; 7) 
of another system with shear rate 7, obeys Eq. (60) at large time (7T >> 1 ), namely 

r Xs t(r,e±;2 7 ) = £U0,£l;7) - D st (r,^ ± ;i)- (61) 

This symmetry of the theory is independent on the model parameters and is therefore 
particularly suited for comparison with other models and for experimental checks. Actually, 
the linearization of the equation of motion provided by the large-A" model is expected to 
provide reliable results above the critical temperature and we expect the relation (61) to be 
observed also in more realistic systems. 

T = T c ( 7 ) 

At T = T c (7) the fluctuation-dissipation plot is the same as for T > T c (j). In the 
following we will focus on the fluctuation-dissipation ratio 

x M= T mm- (62) 

and, in particular, on its limiting value 

Xoo = lim lim X(t,t') (63) 

In systems without drive this quantity takes the value Xoo = 1 for equilibrium systems while 
= in the low temperature phase of coarsening systems with a non-vanishing critical 
temperature [18,21,22]. For critical systems the value of has been computed for the 
random walk and the Gaussian model where the non trivial value = 1/2 is found [23]. 
This quantity has been also computed for some coarsening systems at the lower critical 
dimension d^, such as the X-Y model in d = 2 [23] and the Ising chain [24], where the 
same result Xoo = 1/2 is recovered. However, concerning systems quenched to the critical 
temperature above di, the value of X^ turns out to possess a different value. Numerical 
simulations of the Ising model with heat bath dynamics in d = 2 indicate that X^ ~ 0.26 [25] 
and the exact solution [25] of the spherical model (that is equivalent to the large- iV model) 
gives 

= d —^- 2<d<4 (64) 

Xoo = \ d>4. (65) 
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In [25] it is proposed that this quantity is a novel universal quantity of non-equilibrium 
critical dynamics. 

Under the action of a shear flow, Eqs. (65) is obeyed both for d = 2 and d = 3, as we 
show below. Indeed, concentrating on the case d = 3 for clarity, taking R(t, t') from Eq. (51) 
(with y(t) = const) in the asymptotic limit t — t'^> t' , and evaluating dD/dt' from Eq. (42) 
one obtains 

x °° = l- (66) 

It can be shown that the same result = 1/2 applies also in the case d = 2. Then, with 
respect to the case 7 = where Eqs. (64,65) hold, the effect of the flow is to shift by two 
the effective dimensionality of the system, as already discussed in the previous Section. 



T < T c ( 7 ) 

Summarizing the results of the previous sections, in the regime t w ^> t s one has 

D{t, t w ) = D st {r, oo) + D ag {t/t w ) 



(67) 



Xst(T, OO) 



(68) 



Before discussing the behaviour of the flowing system it is useful to overview the be- 
haviour without shear that is plotted in Fig. 4. In this case a structure like (67,68) is also 
found. 
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FIG. 4. Tx(t,t w ) is plotted against D(t,t w ) for a quench to a final temperature T = T c ( r y)/2, 
with t w = 10 5 , d = 3 and A = 1. Solid lines correspond to 7 = 10~ 2 , 10 , 1, from top to bottom. 
The dashed line is the case with 7 = 0. 
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In the short time separation regime (see Sec. V), or stationary regime, D st {r, oo) decays 
from D st (0, oo) = M 2 - M 2 to D st (r ~ t w , oo) ~ while D ag (t/t w ) ~ £> aff (l) = remains 
constant. Therefore, in this regime D(t,t w ) ranges from Mq down to M 2 . Txstij, oo) grows 
from zero to Txsti^ — t w , oo) ~ Mq — M 2 and is related to D st (r, oo) by the equilibrium 
FDT, Eq. (60). This gives rise to the straight line with negative slope on the right of Fig. 4. 
Notice that this line is nothing else than the fluctuation-dissipation plot of a system in the 
equilibrium state above T c , shifted by the amount M 2 along the horizontal axis. 

In the large time separation sector, or aging regime, D st (r, oo) ~ whereas D ag {t/t w ) 
decays from D ag (l) = M 2 to zero for t/t w large. Hence the aging regime corresponds to the 
region < D(t,t w ) < M 2 of the fluctuation-dissipation plot (Fig. 4). The response function 
stays constant Tx{t, t w ) = Mq — M 2 because the aging part of this quantity is negligible for 
large t w . These behaviours account for the flat part on the left side of Fig. 4. Given this 
structure it is also clear that the whole time dependence of the response can be absorbed 
through D(t,t w ), namely x(Mw) — x(D), as proposed in Ref. [2]. 

Next we turn to the case with shear. In this situation, given Eqs. (67,68) and the 
behaviour of D st {r, oo), D ag (t/t w ) and Xst(j, oo) obtained in Sees. V,VI (Eqs. (42,44,54)), 
most of the considerations discussed for 7 = still apply. In particular in the short time 
separation regime D st (r, 00) decays from D st (0, 00) = Mq — M 2 to zero while D ag (t/t w ) ~ M 2 
keeps staying constant. In this regime D(t, t w ) ranges from Mq down to M 2 and Txstij, 00) 
grows from zero to its limiting value Txstij — t w ,oo). However, when shear is applied, 
the relation between Xst( r ) 00 ) an d D st (T, 00) is not the equilibrium FDT (60) but, instead, 
the relation x(D) of the case T > ^(7) discussed above (shifted by the amount M 2 along 
the D for 7 = 0). Therefore, as discussed already for T > ^(7) in the stationary 

regime, one has both a region r < t s where standard FDT (60) holds (on the right of 
Fig. 4), and a time domain r s < r < t w where the relation x(D) deviates from the linear 
behaviour (60). Finally, for r > t w (namely D < M 2 (7)), the aging regime is explored with 
a flat fluctuation-dissipation plot. 

The pattern of violation of the FDT discussed insofar for coarsening under shear flow 
can be compared with the behaviour of driven mean field models for glassy kinetics [3]. In 
both cases the system without drive exhibits a phase transition at a critical temperature 
T c characterized by an aging dynamics below T c , and a qualitativly similar behaviour of 
the autocorrelation function. However several differences occur between these systems when 
the drive is switched on, some of which are revealed by the fluctuation-dissipation relation. 
The most important difference is due to the fact that mean field glass models always at- 
tain asymptotically a stationary state, regardless of the temperature, while, as discussed 
in Sees. III,V, the model considered in this paper exibits an aging kinetics below a critical 
temperature even in the presence of shear. Furthermore, in glassy models, the fluctuation- 
dissipation relation x(D) is a broken line for every temperature and the two slopes of these 
lines are associated to the existence of two well defined temperatures in the system, the bath 
temperature T and the effective temperature T e ff of the slow degrees of freedom. A simi- 
lar behaviour is found in molecular dynamics simulations of binary Lennard- Jones mixture 
under shear flow [26]. On the other hand, at least in the stationary state observed above 
T c (7), the only one that can be compared with the time translational invariant states of the 
glassy case, the model considered here shows a more complex relation x(D), as it can be 
clearly observed in the inset of Fig. 4. This suggests that, despite some qualitative similari- 
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ties, coarsening systems under the action of an external drive behave quite differently from 
glassy systems and that the violation of the FDT is an efficient tool to detect it. 

VIII. CONCLUSIONS 

In this paper we have studied analytically the out of equilibrium kinetics of a solvable 
model for binary systems subjected to a uniform shear flow. Besides the relevance of this 
subject as a first step for a better comprehension of sheared two-component fluids, which has 
recently been by itself matter of intense theoretical and experimental interest, this model 
is particularly suited for investigating off-equilibrium systems under a general and wide 
perspective. Actually, the properties of non-equilibrium states are an important concern of 
modern statistical mechanics and many efforts have been devoted to widen our knowledges 
in this field. In particular, much interest has been payed to aging systems, namely systems 
which retain memory of the time t w elapsed since they have been brought out of equilibrium, 
such as glasses or spin glasses. The aging properties are usually extracted from the behaviour 
of the autocorrelation function D(t,t w ) which carries an explicit dependence on both times, 
up to the longest observation scales. On the other hand, the so called driven systems, 
to which external energy is pumped from the outside, generally lose memory of the initial 
condition and, after a microscopic time, a time translational invariant state is entered, which, 
because of the external drive, is not an equilibrium one, in the Gibbs sense. It has been 
argued that an important piece of information on non-equilibrium states, both aging and 
stationary, is encoded into the linear response function to an external perturbation, and, in 
particular, in the so called off-equilibrium generalization of the FDT, the relation between 
the response function and D(t,t w ). Hopefully, this non-equilibrium linear response theory, 
should bear the amount of relevant implications and the generality of the equilibrium case. 
However, nowadays, our understanding of this important issue is still incomplete and further 
investigations are necessary in order to deepen our insight. 

The model we have studied in this paper is particularly suited for this analysis because 
it offers both a framework where exact calculations can be carried over and a rich phe- 
nomenology due to the presence of both stationary and aging off-equilibrium states. These 
two behaviours are separated by a phase transition at a critical temperature T c (t) which 
increases with 7 and gives a transition also in d = 2. We stress that this is a transition be- 
tween non-equilibrium states. Nevertheless, in this model the existence of a non-equilibrium 
parameter 7 allows one to compare the properties of the transition with the well known 
equilibrium case with 7 = 0. In doing that, one discovers that a wealth of similarities exist 
greatly helping the analysis of what goes on out of equilibrium. In particular, in the low 
temperature phase, where the aging behaviour coexists with the drive, the autocorrelation 
function can be splitted into a stationary and an aging part, along the same lines as for 
7 = 0. The former describes the correlation between fast degrees of freedom while the lat- 
ter, which decays as a power law with a generalized Fisher-Huse exponent, is due to the 
slow variables responsible for the aging properties. This analysis allows a better compre- 
hension of the violation of the FDT and the recognition that a mechanism very similar to 
that observed without shear is at work. This provides the basis for understanding why the 
fluctuation-dissipation plot of Fig. 4, with the flat part characteristic of corsening systems, 
comes about also in sheared systems. 
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In aging systems without drive a definite progress in generalizing the FDT out of equi- 
librium has been achieved by realizing that certain fundamental properties of the system are 
encoded into the fluctuation-dissipation relation. A first step in this direction was made by 
interpreting T e ff = TX _1 (C) as an effective temperature of the system, different from that 
of the thermal bath. Furthermore, it was also shown by Franz, Mezard, Parisi and Peliti [27] 
that the relation between response and autocorrelation function bears informations on the 
properties of the equilibrium state towards which the system is evolving. These advances 
qualify the fluctuation-dissipation relation as a fundamental one also far from equilibrium, 
providing an important tool in this difficult field. However, a generalization of these concepts 
to off-equilibrium driven systems is lacking, particularly for the low temperature phase of 
the model considered here, where aging coexists with the drive. Further studies are therefore 
needed to discover if and which properties of the non-equilibrium state are encrypted into 
their fluctuation-dissipation relation. 



IX. APPENDIX A: THE LAPLACE TRANSFORM OF F(T) 

In this Appendix we compute the Laplace transform of f(t) defined as 

f L (s) = / e- ts f(t)dt (69) 
J o 

where fit) is given in Eq. (32). In doing so we will obtain also an explicit form for the 
quantity P(x) defined in Eq. (18) due to the relation 

P(x) = 2f L (2x). (70) 

This relation can be easily derived by inserting the expression (31) into Eq. (69) and com- 
paring the result with the definition (18). 



d=2 

We first calculate the Laplace transform of the function fit) of Eq. (31) in d = 2. It can 
be written as 

fHs) = r /(*>* = - f ( 4 - 7 ^ + "7^1 "* dy (71) 

J K ' Jo JK J 4tt Jo 1 + Ay \ (1 + Ay) 2 3(1 + Ay) J y K J 

e ^A yoo e -z / [z^A/s - 1)V 4 (z-fA/s - l) 3 s\~' 2 
dz 4 — - - -|- — - 



4tt z \ A^ 2 z 2 3 A^z 



where the variable z = s(l + Ay)/( r yA) is introduced and A is defined in Eq. (21). The last 
integral can be split into two pieces: The first is given by 



87r J ix Z x /l + h(z) 

where the function 



/ dz- r —= (72) 
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{z^A/s-lfs 2 f Z1 A 1\ 

can be shown to be less than one in the integration interval. Then the square root of Eq. (72) 
can be expanded as a power series and integrated term by term. The result is 



with r n = (— l) n (2n — l)!!/(2n)!!. The other contribution to f L (s), neglecting all terms of 
order 0(s 2 , A" 1 ), is given by 



oo „ „ oo 



p ~«A poo p * ; 

e dz - (75) 

7 V 1 + 12? 

Proceeding along the same lines as for the previous integral (72) and adding the result to 
expression (74) one finds 

1 / „ oo -1971+I/2 2n+l (In -I \ A /T9o\ 

/'(.) = I (NVT2A) + ! + !>„ + + E li_f^_ (n ^ .n±fi) 

+ 0(s 2 ,A- 1 ) (76) 

where Do = E^i^^Si) = ln2 ~ 1 ~ -°- 3068 > D 1 = C - 2 - D ^ -1.7296 with 
C ks 0.577216 being the Euler's constant. Eq. (76) can be written more compactly as 

M 2 

/L(S) ^ 27^) -^ + 5 2 alo ga (77) 
where only the first term in the series of Eq. (76) is retained and 

M 2 1 

/*(0) = ^ = ^/l2A) + ln2) (78) 



A 2 = -^V(£>! - ln()/l27^)) ;fi 2 = ^ 7 - 1 - (79) 



d=3 

Similarly to the case d = 2, starting from Eq. (32), f L (s) can be written as 



f L, \ _ V^ A3/2 r <Q1 (a _ y' A2 4 y 3A \ 5 w , 8m 

/3 W ~ 4(271)3/2 io (l + Ay)^^ 4 (l + A 2/ ) 2 + 3(l + Ay) y ' ^ l8UJ 

r [°° j e ~ Z (a (z-fA/s - 1) 4 s 2 A(z^A/s-lfs\~^ 2 

S I CLZ -r-rr 4 — - - -|- — - 



4(2tt) 3 / 2 v z 3 / 2 V A^ 2 z 2 3 A 3 1Z 
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The last expression can be again calculated splitting the integration interval. A first contri- 
bution comes from the term 

y/s fVTZsh e -z I 

dz- 



2 z 3 / 2 



V 1 + h{z) 



Adding this contribution to the other term 



VI f°° , e~ z 1 >fi fj, 

2 J^* ,?3/ 2 /. . ^2 0/10U/4 2j 



7 



V^Tfg 2(12)V4^ 2n + 3/2 



oo p/i /n\ oo /2n+l o \ ion+1/2 

-^S^^Sfn^)^— + ^,,-v, (82) 

one obtains 

1 r /— r x/7 i/4 s r(l/2) ~ Z 2 ^ 1 2 \ 2n , , /2 12 n+1 / 2 



/ L (s)= 



4(2tt) 3 / 2 



12 V4 ^ 1 2 ^oVto 2fc + V 7 2 " +1 



+C(s 2 ,A- 1/2 ) (83) 

where T = S^r W ( ^+ff_ 1/2) « -0.847, 7\ = £~ 0^5^175 « L854L The behaviour 
for small s is the only relevant for the asymptotic dynamics, as discussed below Eq. (93). 
This is given by 

/L(s) " mn) ~ A3S + 5358/2 (84) 



with 



/i( °> = 2^) = i(2^<V^+ I^T^o) (85) 



A 12V4 T l D V^12 _! , , 

^"4(2^Vf 5 3 = W^ 7 • (86) 



It can be shown that the coefficients To, T\ can be also written in terms of the T-function as 

T _ r 2 ( 3 /4) T _ r 2 (i/4) 



X. APPENDIX B : THE FUNCTION Y(T) 

We want to compute y(t), the inverse Laplace Transformation of y L (s) of Eq. (35), 
defined by 

y(t) = — ds y L (s)e st (87) 

Z7T1 J a— too 
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where a has to be chosen in such a way that all the poles of y L (s) are on the left of the path 
of integration. 

We first study the pole structure of y L (s) in the half-plane x > with s = x + iy. The 
real part R(x,y) of the denominator of y L (s) is given by 

R(x,y) = x -2r -4TRe{f L {x,y)) (88) 

where 

Re (f L (x, y)) = J™ dtf(t)e- xt cos(yt). (89) 

Given that f(t) > for any time, (see Eq. (32)), for non-negative values of x, Re(f L (x ) y)) 
is a monotone decreasing function of both x and y. This implies, from Eq. (88), that R(x, y) 
is a monotone increasing function in the half-plane considered. As a consequence, R(x, y) 
remains positive for all x > if R(0, 0) > 0. Since from the definition of T c (pj) given in 
Eq. (23) one has 

i?(0,0)>0, T<T C ( 7 ) (90) 



i?(0,0) = 0, T = T C ( 7 ) (91) 

i?(0,0)<0, T>T C ( 7 ), (92) 

for T < T c (7), R(x,y) is always positive for x > and y L (s) has no poles in the complex 
positive half-plane. In particular at T = T c ( 7 ) there is a pole in the origin. On the other 
hand, for T > T c ( 7 ), the monotonic behaviour of R(x, y) and Eq. (92) imply the existence of 
a positive p such that R(x, y) > R(p, 0) = for any x > p. Since the imaginary part of the 
denominator of y L (s) is zero everywhere on the real axis, one can conclude that the largest 
pole is located on the real axis at p and, in the Bromwich contour (the vertical integration 
path of Fig. 5), a > p. 

We can now calculate explicitly the integral of Eq. (87) in the different temperature 
regimes. 



T<T C ( 7 ) 

In this case there are no poles in the positive complex plane and the Bromwich contour 
can coincide with the imaginary axis. The theorem of residues can be applied choosing a 
closed contour in the negative complex half-plane. Due to the non analytical part of f L {s), 
there is a branch cut in the complex plane along the negative real axis. Then we choose 
the contour of Fig. 5 and the only non vanishing contributions come from the Bromwich 
contour itself and from the integral along both sides of the branch cut. Then, one has 

J dsy L (s)e st = 2my(t) + J™ dxe~ xt (y L (s) \ s=xei « -y L (s) \ s=xe -^) . (93) 

Since we are interested in the large time behaviour of y(t), we can compute the above 
integrals taking only the small s behaviour of f L (s) given in Eqs. (77,84). We obtain 
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jf dsy\s)e st = 2iny(t) - 2ttz^| (t + A M 2 ) t" 2 , d = 2 (94) 

| e dsy L (s)e st = 2my{t) - ( T + A ° M o) ^ d = 3 ( 95 ) 

The above expressions are equal to the sum of the residues of all eventual poles inside G. In 
the negative complex plane these poles can only give contributions exponentially decreasing 
with time which are asymptotically negligible with respect to the power law behaviour of 
Eqs. (94) and (95). Henceforth we obtain expressions (36,37) in d = 2 and d = 3 respectively 

▲ 

Im(s) 




— ► 

Re(s) 



FIG. 5. The closed contour adopted to compute the integral (87). 



T = T C ( 7 ) 

The same contour of the T < T c (jy) case can be chosen since the pole located in 
the origin is outside 0. The evaluation of the integrals proceeds differently in d — 2 and 
d = 3. In d = 2 the computation strictly follows the one presented for T < T c (j). The 
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non vanishing contributions on still come from the Bromwich contour and the integrals 
around the cut. One obtains 



L ds ^ a)e " - - 2 "2ftk> 6 + ti) (96) 

which gives, at the end, the result of Eq. (36). In d = 3 the contribution coming from the 
small half-circle around the origin gives a constant term. Since, as before, the integral along 
both sides of the cut has a decreasing power law behaviour, asymptotically, the contribution 
from the half-circle around the origin is the dominant one and 

/ 9 ^ We ^2„ W -„i±if^M (97) 

resulting in Eq. (37). 
T > T c ( 7 ) 

In the case of a temperature T = Tj^f) + 5T with 8T/T C {^) <C 1 the preliminary 
discussion about R{x, y) assures the presence of a pole on the positive real axis close to the 
origin. The location (p, 0) of this pole is given by the solution of the equation 



R(p, 0) = (98) 
that, for p -C 1, as 5T/T C (^) <C 1 implies, reads 

p\og(p) = M Q 2 5T/(T c ( 1 )2B 2 T), d = 2 (99) 

p(l/2 + 2TA 3 ) = M*5T/T c ( 1 ), d = 3. (100) 



At this point a similar procedure as in the case T = T c (j) can be followed giving the 
asymptotic behaviour for y(t) of Eqs. (36) and (37) in d — 2 and d — 3, respectively. 

XI. APPENDIX C: ASYMPTOTIC BEHAVIOUR OF THE ZERO MODE OF 

C(K,T) FOR T < Tcin) 

The value C(k = 0,t ) can be obtained from Eq. (28) and reads 

An 2T /•*« 

C(k = 0,t ) = ^ + — y(z)dz. (101) 

y{h> y{to) Jo 

The integral appearing in the above equation can be splitted into two terms: 

rto fco poo 

/ y{z)dz = / y(z)dz - / y(z)dz. (102) 

J0 J0 Jt 

From Eq. (35) one has 
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1 ]\yf2 A 

y(z)dz = y L (0) = + ^^-r (103) 
yK J y K ' 2M 2 M 2 2T C ( 7 ) v ; 

where we have used Eqs. (78,85), while the asymptotic behaviour of y(t) (Eqs. (36,37)) gives 

r y{z)iz = _JL_ VUT+^MS 2 dP 

Jto yV ; (8vr) d / 2 7 M 4 d V ; 

Putting together Eqs. (101,103,104) one finally obtains 



r(h c\t\ 1 ( M 2 4T V12 .- d/2 \ T + AoK 

C(A; = °' to) = { ~ WW 2 V° J • (105) 

Notice that the dependence on the initial condition, namely the parameter A , drops out 
for t 3> t s in Eq. (105), by virtue of the expressions (36,37) of y(t ), because the scaling 
regime is entered and the memory of the intial state is lost. 



XII. APPENDIX D: THE AUTOCORRELATION FUNCTION 

The evolution equation (40) for the autocorrelation function D(t,t') can be integrated 
starting from a generic time to chosen arbitrarely between the time of the quench (t = 0) 
and the observation time t. The formal solution for D(t,t') then reads 

D(t,t') = D T (t,t') + D (t,t') (106) 

where 

DtM = - 2T / ^L e -(^(7)/2)/A2 f* dze -2jr>CHu) d u e jr'^^y(z) (107) 

y/y(t')y(t) J ( 27r ) Jt ° 

and 

r» (*,0 = , y{t ° ] I ^C*(^-t ),t )e-( fc2+d W^^ (108) 

with d( 7 ) = K 2 (t-t')-k 2 . 

We will take the time t large enough (t 3> Ts) that y(t) is given by Eqs. (36,37). 
Furthermore it can be shown that, due to the condition (21), the function d{^) can be 
neglected in Eqs. (108,107). 

We will first consider the function £)y(t, t'). Carrying out the integral one obtains 

AT 

D T (t,t') = x (109) 

/ y(z)(t m - z + r M y d/2 ^= ^dz 

Jto . H 2 ((t-z)2- T 2/ 2 )2 4 2 (( t -^)3_ r 3 /2 ) 

' (t m -Z+T M ) 2 ' 3 ' (tm~Z+T M ) 



27 



with t — t — t' and t m — (t + t')/2. The role of the cutoff A is crucial in the factor 
(t rn — z + tm)^ 4 ^ 2 in order to avoid divergences at z = t m ; however, for large t m > r s , one 
finds that it can be neglected in the two terms under the square root. Then, introducing 
the variable u — 1 — z/t m , one has 



D T (t,t') 



4T t~ d/2 fl-to/tr, 



where 



p(u, r, t m ) 



(8vr) d / 2 7 Jr/(2t m ) 



p(u,r, t m )Q(u,T, t m )du 



(110) 



(111) 



and 



0(u,T,t m ) = < 



f e -2^t mU ; for T > T c ( 7 ) 

1 ; for T = T C ( 7 ), d= 3 

Vlog(t +r/2)lo g(f for = = 

log[t m (l-u)] _ c\ ;/) 

5 " : lor V < T r (" ) . 



(112) 



I \y/l-{T/2tm) 2 



Turning now to the function D (t,t'), in the limit t'/t 3> 1, due to the presence of the 
exponential factors in Eq. (108), Do(t,t f ) can be evaluated by approximating the one-time 
correlation function C (JC{t — t ),t ) at the time to entering Eq. (107), whose expression is 
given in Eq. (28), with its value at k — 0. Then one has 



D (t,t') = 

C(k = 0,t )y(t ) {t _ r/2 + 1/2Ard/2 



(113) 



„ 2 (t 2 -r 2 /2)2 4 2 (t 3 ~r 3 /2) 



/4 _ ~2 ^ -■ /jV |_ ±^2 

" • ' (t-r/2+l/(2A 2 )) 2 ^ 3 ' (t-r/2+l/2A 2 ) 

From this point on one must distinguish different temperature ranges. 



T > T c ( 7 ) 

Starting with D T (t,t'), from Eqs.(110,lll,112), one has 



D 



2\-d/2_ 



4 + § 7 V + t r* (l - g) 



dy. (114) 



In the limit t m — > oo, D T (t m , r) becomes a TTI quantity 



a, 



4 + § 7 V + ^ r* (l - 



dy. (115) 



With respect to D (t,t'), by inserting the asymptotic behaviours of y(t) ( Eqs. (36,37)) 
in Eq. (113), one immediately obtains 
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D (t,t') ~ e - 2 ^ t+t '^ 2 (116) 

From expressions (115,116) it is clear that, in the limit t' ^> D (t,t') is negligible 
with respect to D T (t,t') so that 

D(t,t)c«D H (T,t ± ) (117) 
and D st (r,£±) is given in Eq. (115). 



T = T c ( 7 ) 

In this case the same results (114,115) are found for D T (t,t') but with ^_ L 1 = 0. In 
particular the large t m behaviour of D(t,t') is given by 

AT r°° , 1 

D st (r,oc) = -— ^ (y + l/2AY d/2 , — dy (118) 



4 + § 7 V + (l - 



with corrections of order t m d A 

Regarding D (t,t'), from Eqs. (36,37) one has 

A)(M')~(^) (119) 

with logarithmic corrections in d = 2. 

Comparing Eq. (118) with Eq. (119) one can show that again, in the limit t' 3> t s , 
D (t,t') is negligible with respect to D T (t,t'), as for T > T c ( 7 ). Therefore 

D{t,t') ~D st (r, oo). (120) 

The expression for D st (r, oo) is Eq. (118). 



T < T c ( 7 ) 

For T < T c (7) we study the behaviour of Drp(t,t') separately in the temporal regimes 
r/t m <C 1 (quasi-stationary regime) and r/t m ^> 2/3 (r/t' ^> 1) (aging regime). In the 
quasi- stationary regime, the function p(u,r,t m ) of Eq. (Ill) is of order t^ +2 ^ 2 for u <C 1 
and decreases to zero as p(u, r, t m ) ~ u ~^ d+2 ^ 2 when it > l/ 7 t m . In the same regime, 6 ~ 1 
when m C 1 while it diverges like (1 — w)~( d + 2 )/ 2 for u — > 1 — t /t m - In the limits 7 t m — > oo 
and to Am. - *■ 0, the leading contributions to the integral of Eq. (110) coming from the 
singularities of the functions p(u, r, t m ) for u — > and of 0(it, r, t m ) for it — > 1 — to Am give 

J Dr(t m ,r)= (121) 



4T t^ 2 

(87T) d / 2 7 

D ^)+(^" /2+ °fc i/2 >- 



/■oo rl~t /(2tm) 

/ dup(u,T,t m ) + V3 / duQ(u,T,t m ) 

JrlVltm) JO 



+ 0(C /2 ) 
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Similar considerations applied to the aging regime imply 
D T {t m ,r)= (122) 

4T VT2 t _ d/2 ,ty<? + y 2^ + 



(87r) d / 2 c?7 VW l A 2-(i-t'/t)3 o (2-(i-tyt) 2 ) 2 

Y l+t'/t J (i+t'A) 2 

For D (t,t') one finds 

A)(M')= (123) 
/ 2 4T >/l2 _ d/2 \ (t\-*P„ 2^ 



(87r) d / 2 d-f J \t'J V 7 ; " /, 2-(l-f /t)3 (2-(l-t'/t) 2 ) 2 

Y i+*7* (i+t'A) 2 

— # 

where the expression of C(k — 0, to) obtained in appendix C has been used and the limit 
t' ^> r s is assumed. 

Summarizing, summing up D (t,t') (Eq. (123)) with D T (t,t f ) (Eq. (121)), in the quasi- 
stationary regime one obtains 

D(t,t')~D st (T,oo)+M 2 (124) 
with D st {r, oo) given in Eq. (118) while, in the aging regime, from Eqs. (123,122) one has 



V ; 9W ; UV V ' ' / /1 2-(i- f yt)3 (2-q-tyt) 2 ) 2 

V 1+f/t ° (i+t'/t) 2 

thus defining the quantity D ag (t/t'). 
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